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REPEATED  MEASUREMENTS  DESIGNS,  III 
BALANCED  CONNECTED  DESIGNS 


By 

C.  Magda  and  A.  Hedayat 

Department  of  Mathematics 
University  of  Illinois,  Chicago 

ABSTRACT 

Repeated  measurements  designs  (RM  designs)  are  experiments 
in  which  an  experimental  unit  is  exposed  to  a sequence  of 
treatments  during  a number  of  periods  (days).  A measurement 
taken  on  a unit  depends  on  the  treatment  administered  in  the 
previous  period  (residual  effect),  the  treatment  given  during 
the  present  period  (direct  effect),  the  unit  itself  (unit 
effect)  and  the  present  period  (period  effect). 

In  the  first  part  of  the  paper,  after  defining  the  con- 
cept of  balanced  and  connected  RM  design,  we  proceed  in  cons- 
tructing large  families  of  balanced  and  connected  RM  designs 
when  the  number  of  periods  is  less  than  the  number  of  treat- 
ments. In  the  last  section  we  discuss  the  structure  and  the 
ranks  of  the  C-matrices  for  direct  and  residual  effects  which 
arise  from  these  designs,  as  well  as  the  form  of  the  best 
linear  unbiased  estimators  of  estimable  contrasts  of  residual 


and  direct  effects 
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REPEATED  MEASUREMENTS  DESIGNS,  III 
BALANCED  CONNECTED  DESIGNS 


C.  Magda  and  A.  Hedayat 

Department  of  Mathematics 
University  of  Illinois,  Chicago 


1 . Introduction 

A statistical  design  in  which  experimental  units  are  used 
repeatedly  by  exposing  them  to  a sequence  of  treatments  is 
called  a repeated  measurements  design  (RM  design).  In  a RM 
design  we  administer  t treatments  during  p periods  to  n 
experimental  units . The  parameters  t and  p are  usually 
known  to  the  experimenter.  RM  designs  have  been  considered 
in  the  literature  under  a variety  of  names  such  as  crossover 
or  changeover  designs,  time  series  designs  or  before-after 
designs.  In  these  designs,  apart  from  the  direct  treatment 
effects,  residual  or  "carry-over"  effects  are  usually  assumed. 

It  has  been  shown  by  A.  Hedayat  and  K.  Afsarinejad  (1978), 
that  if  certain  balance  for  the  residual  and  direct  effects 
is  assumed,  these  designs  become  optimal  in  a large  class  of 
designs.  The  cyclic  RM  designs  that  we  shall  construct 
(precise  definitions  will  follow),  apart  from  having  these 
balancing  properties,  are  also  connected,  i.e.,  all  contrasts 
of  residual  and  direct  effects  are  estimable.  Their  cyclic  struc” 
ture  makes  them  easy  to  implement  in  practice.  We  also  give 
the  best  linear  unbiased  estimators  (b.l.u.e.'s)  for  some  of 
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the  contrasts . Tables  of  such  designs  for  small  values  of 
t and  p can  be  found  at  the  end  of  the  paper.  The  authors 
hope  that  they  would  suffice  for  most  practical  applications. 


2 . Definitions  and  Model 

Let  us  present  the  set-up  in  which  we  will  be  working. 
Suppose  d is  a two  dimensional  array.  If  a is  in  position 
d(i-l,j)  in  d and  b is  in  position  d(i,j)  in  d we 
say  that  a precedes  b (in  the  order  of  application). 

We  say  that  a follows  b if  b precedes  a. 

Definition  2.1.  A balanced  RM(t,p,n,X^,X2)  design  is  a 
pxn  array  d whose  entries  involve  t treatments  (labeled 
0,1, . . . ,t-l) , satisfying  the  following  conditions: 

(a)  the  treatments  along  any  column  are  distinct. 

(b)  A treatment  occurs  exactly  X^  times  in  each 
row  of  d . 

(c)  In  the  order  of  application,  each  treatment  is 
preceded  in  d by  each  other  treatment  exactly 
Xg  times . 


Part  (a)  in  the  above  definition  has  been  included  mostly 
because  of  esthetic  statistical  reasons,  such  as  treating  a 
column  in  the  array  as  a fixed  size  block,  etc...  We  refer 
to  the  treatment  in  the  (i,j)th  entry  of  the  design  d 
by  d(i,j).  The  basic  divisibility  conditions  for  the  existence 
of  a balanced  RM(t,p,n,X1,X2)  design  are  the  following: 


rr 


n = X^t  and  X1(p-1)  = X2(t-l). 

Definition  2.2.  a shift  starting  at  i of  (0, 1,2 , . . . ,t-l) 
is  (i,i+l, . . . ,t-l,0,l, . . . ,i-l) . 

Definition  2.3.  We  say  that  two  shifts  of  (0,1,2, .. .,t-l) 
are  distinct  if  they  start  with  different  symbols. 

Definition  2.4.  A balanced  RMftjPjnA^Xg)  design  is  called 
cyclic  if  it  can  be  decomposed  into  X-^  subarrays  of  size 
Pxt  and  the  rows  of  each  pxt  subarray  are  shifts  of 
(0,1,2, . • • , t— 1 ) • 

Throughout  the  paper  the  underlying  model  will  be: 

Ylj  - » + “i  + + Td(l,J)  + »cl(l-l,j)  + 

where  e^j  are  assumed  to  have  mean  zero  and  covariance  matrix 
a2I.  The  overall  effect  is  p,  is  the  effect  of  the  ith 
period,  is  the  effect  of  the  Jth  experimental  unit, 

Td(i  j)  is  the  direct  effect  of  treatment  d(i,j)  and 

pd(K  -t)  is  tile  residua^  effect  of  treatment  d(K,t).  During 
the  first  period  (l.e.,  for  i=l),  no  residual  effects  are 
assumed,  so  j)  are  zero  for  J* 


4. 


3.  Construction  of  Balanced  hlil  Designs 
with  minimal  number  of  units 


In  this  section  we  construct  families  of  cyclic  balanced 
RM  designs.  The  main  result  is  contained  in  Theorem  3.1. 

We  begin  with  a couple  of  useful  propositions.  As  we  shall 
see  in  the  proposition  that  follows,  the  advantage  of  con- 
structing balanced,  RM  designs  by  shifts  is  that  condition 
(c)  of  Definition  2.1  needs  to  be  satisfied  for  one  treatment 
only.  Let  the  treatments  be  {0, 1,2, . . . ,t-l J . 

Proposition  3 »1 . A pxn  array  d,  with  X ^ subarrays  of 
size  pxt  each,  is  a cyclic  balanced  RM(t,p,n,X1,X2)  de- 
sign if  all  the  rows  in  each  pxt  subarray  are  distinct 
shifts  of  (0,1,2, ... ,t-l)  and,  in  d,  the  treatment  0 is 
followed  in  order  of  application  by  every  other  treatment 
exactly  X 2 times . 

Proof:  Since  the  shifts  in  each  pxt  subarray  are  assumed 
distinct,  (a)  is  satisfied.  Clearly  (b)  holds.  To  prove  (c), 
fix  treatment  j.  By  assumption,  0 precedes  j exactly  X2 
times.  So  let  i be  different  of  0 and  j.  The  treatment 
j-i  (mod  t)  is  nonzero  and  hence,  preceded  by  0 X2  times 
in  d-  There  are  therefore  exactly  X2  rows  of  length  t 
in  various  pxt  subarrays  in  which  j-i  (mod  t)  is  preceded 
by  0.  In  all  these  rows,  J is  preceded  by  i,  because  each 
such  row  is  assumed  to  be  a shift  of  (0,1,2, ... ,t-l) . There- 

L Zl 
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fore  j is  preceded  by  i exactly  Xg  times. 

The  following  proposition  reduces  the  construction  of  a 
cyclic  balanced  RM  design  to  a partition  problem  and  is  our 
tool  for  constructing  these  designs. 


Proposition  3 .2 . Let  A.^p-1)  = (t-l)  (i.e.,  Xg  = 1). 

A cyclic  balanced  RM(t ,p,n,X l)  design  d exists  if  and 

X1 

only  if  {1,2,3, ...,t-lj  = U P^  (disjoint),  |Pj  = p-1 

i=l  1 1 

such  that  for  some  ordering  of  the  elements  in  Pj^ , the 
successive  partial  sums  of  these  elements  are  all  distinct 
and  nonzero  ( modulo  t ) , for  all  i,  1 < i < \ 


Proof:  For  simplicity  of  notation,  in  the  proof  that  follows 
we  will  not  use  the  d(i,j)  notation  introduced  earlier. 


Note  that  the  condition  on  the  successive  partial  sums  in  P.^ 

to  be  distinct  and  nonzero  (mod  t)  can  be  expressed  equlva- 

K 

lently,  by  subtracting  two  partial  sums,  as  E s.  4 0 (mod  t) 

t=j 

for  all  1 < j < k < p-1,  where  P.^  = {(s.j,)j. 


Suppose  we  have  a cyclic  balanced  RM(t,p,n,X^,l)  design 

4*  h 

d.  Let  f^j  map  the  entries  of  the  row  into  the  corres- 

1.  u i.  v- 

ponding  entries  of  the  (j+l)  row,  in  the  i1'  pxt  subarray 
of  d.  Because  d is  a cyclic  design,  f^j  has  the  following 


property:  fij(x)  = + x*  Let  sij  = Define 

Pi  = ^sij^l<j<(p-l)  = f(sii»s12,Si5,...,Si(p_i))i. 

X1 

U P.  = [1,2,3, • • .,t-lj  holds  because  0 is  preceded  in  d 
i=l  1 

by  everything  else  exactly  once.  On  the  other  hand,  in  the 

i.  u 

i^  pxt  subarray  of  d,  the  entries  along  any  column  are 
distinct,  so  in  particular, 

f’ll(D) » f j_2  ( f j_i(  o) ) > • • Pi(  p-2 ) ( • • • ( ))•••))  are 

distinct.  But  j_i)  ( • • •(fn(°))  • • • ) ) = sil+si2+*  * *+sij 

k 

So  E s,,  4 0 (mod  t)  follows  (by  subtracting  two  such 
t=J  1 

partial  sums),  for  all  1 < J < k < (p-l). 

Conversely,  assume  Pi  = { (sii'si2»si3»  • • • 'si(p-i)  )i  • 
Construct  the  cyclic  balanced  RM(t,pJn,X1,l)  design  d as 

X L- 

follows : In  the  1°  pxt  subarray  let  the  zeros  be  followed 
by  the  sij's  in  the  same  order  as  they  appear  in  P^.  Then 
complete  the  subarray  by  shifts.  Down  a column  (without  loss 
the  first  column)  we  have  0,  8^,Sj+s2,  • . .,s-j+s2+. . .+s^p_1j 
which  are  all  distinct  by  assumption.  Since  we  complete  the 
subarray  by  shifts,  the  entries  along  any  column  of  the  sub- 
array will  also  be  distinct.  With  the  help  of  Proposition  3.1 
we  conclude  the  proof. 
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Remark.  An  analogous  result  to  Proposition  3.?  has  been  estab- 
lished in  which  the  restriction  = 1 is  not  imposed.  The 
authors  chose  to  include  the  present  version  because  of  its 
direct  impact  on  the  results  that  follow  and  because  of  its 
simplicity. 

We  introduce  two  examples  to  illustrate  how  a partition* 
as  described  in  Proposition  3-2,  can  be  associated  with  a 
cyclic  balanced  RM  design. 


Example  1.  Let  t=9*  P=3*  n=36,  X-j=4  and  X2=l.  Consider 
P-L  = ( (1*7)  J*  P2  = C (2,6)  i*  P5  = ( (3*5)  i and  P4  = [(4,8)j. 
These  P-^'s  satisfy  the  requirements  of  the  previous  proposi- 
tion. The  associated  cyclic  balanced  RM(9*3*36*4,l)  design  is 


,012345678 
^ 2 3 4 5 6 7 8 ,0 

80123456  V7 

/0  12345678 

I3  4 5 6 7 8 ,0  1 2 

8012  34  \5  67 


,0  12345678 

U 345678  /0  1 
8012345  '6  7 

,0  12345678 

U 56  78  ,ol2  3 
34567  ^8  012  . 


Example  2.  For  t = 9,  P = 5,  Xx  = 2,  let  P^^  = {(l,7,3,5)j 
and  P2  = ((8,2,6,4)j.  This  partition  corresponds  to  the 
design: 
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,0  1 2 3 4 5 6 7 8 ,0  1 2 3 4 5 6 7 8 

'a  2345678, 0 ^8, 01234567 

8, 0 123456  W 1 '2  345678, 0 

2 4 5 6 7 8 ,0  1 7 8 ,0  1 2 3 4 5 U 

7801234  V5  6 23  \4  567801. 

Note,  however,  that  the  partition  P1  = ((1, 7,3,4)]  and 
P2  = {(2, 6,4,8)]  does  not  generate  a cyclic  balanced  RM 
design  in  the  same  fashion.  This  happens  because,  in  Pg, 

6 + 4 + 8 = 0 ( mod  9 ) • 

In  the  next  theorem  we  give  a constructive  proof  of  the 
existence  of  cyclic  balanced  RMCt^n,^,!)  designs  for  all 
sets  of  parameters  satisfying  p < t and  the  basic  divisibility 
conditions.  It  will  be  shown  later  that  all  these  designs  are 
connected  (i.e.,  all  linear  contrasts  of  direbt  and  residual 
effects  are  estimable). 

Theorem  3»1«  The  necessary  and  sufficient  conditions  for  the 
existence  of  a cyclic  balanced  RM(t ,p,n,X^, l)  design  with 
p < t are  n = X^t  and  X^p-l)  = (t-l). 

Proof ; (by  construction). 

We  first  give  the  partitions  that  we  use  to  construct 


the  designs,  for  all  cases. 





9* 

case  1.  (p  even),  consider  the  following  partition  of 
* t-l  J : 

Pi  = ((i,(t-l)-ri,  i+Xr  (t-l)-(i+X1),  i+2X1,  (t-l)-(L+2\1),.. 

i+j\r  (t-l)-(i+JX1),...,i  + tejp-  X1)  J for 
0 < J < (pg2),  1 < i < Xx,  and  for  i = X1 
Px  = { (Xx,  (t-l)-X1,  2X1#  {t-l)-2\lf  3XX,  (t-l)-3X1>... 

• - . , JX  -jl,  (t-l)-jXx, ...» (t-l) ) J where  1 < J < (Pg2)  . 

case  2 .1.  (p  odd  and  X1  > 2).  Define  the  partition  as 
follows : 

P±  = { (1,  (t-l)-i,  i+X^  (t-l)-(i+X1),  i+2X1#  (t-l)-(i+2X1),.. 

i+jx^  (t-l)-(i+JX1),...,i  + i^2)x1,  (t-l)-(i  + i£^2lxi 
for  0 < j < ^-Pg^  > 1 < 1 < and  for  i = 

pXi  = {(X1,  (t-l)-Xr  2Xr  (t-l)-2X1,  3XX,  (t-l)-3X1,  ... 
...,JX1,  (t-l)-JX1# ... ,l£^ilxi,  (t-l) ) j for  1 < j < 

case  2 .2 . ( p odd  and  X = 2 ) . 

Let  Px  = {(Mt-l)-l,  3,  (t-l)-3»  5,  (t-l)-5, . . . 

...,  1+2 j , (t-l)-(l+2j),...,l  + 2i£§2l  (t-l)-(l  + 2-^pi))  J 

for  0 < J < and 
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P2  = {(t-1),  2,  (t-l)-2,  4,  (t-l)-4,  5,  (t-l)-6, . . . 
...»  2J,  (t-l-2J,...,(p-l)) j for  1 < j < 
and  (p-1)  = . 


The  proof  consists  in  showing  that  the  conditions  stated  in 
Proposition  3.2  are  satisfied  for  all  the  partitions  above. 

In  all  these  cases  it  is  straightforward  to  check  that 

*1 

U Pi  = (1*2,3, . . . ,t-lj  (disjoint)  and  that  |P.|  = p-1 
1=1  1 

for  all  1 < i < 

k 

To  check  the  condition  E s,  4 0 (mod  t)  for  all 

1= j 14 

1 < J < k < (p-1),  where  P±  = ((s^)^  < ^ < (p_x)  is  the 

Pj^  defined  above,  we  note  that  any  two  (suitable)  neighbor- 
ing entries  in  P.^  always  sum  up  to  t-1,  for  all  1 < i < X^ 
Th*s  fact  enables  us  to  show  that  an  arbitrary  sum  E s.  , 

as  above,  is  never  a multiple  of  t.  A detailed  verification 
for  all  cases  is  given  in  the  Appendix. 

In  practice,  the  number  of  treatments,  t,  and  (possibly) 
the  number  of  periods,  p,  are  known  to  the  experimenter.  It 
is  usually  of  interest  to  construct  a balanced  RM  design  d, 
with  known  t and  p,  in  which  the  number  of  units,  n,  is 
minimal.  This  can  be  achieved  (as  the  divisibility  conditions 
show),  if  and  only  if  X^  and  X2  are  relatively  prime.  The 
designs  in  Theorem  3.1  have  minimal  n since  X2  = 1. 
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It  is  easy  to  extend  the  validity  of  the  previous  theorem 
in  the  following  sense: 

Corollary  3.1.  The  necessary  and  sufficient  conditions  for  the 
existence  of  a cyclic  balanced  RM(t,p,n,X^,X2)  design  with 
p < t and  Xg  |X i are  n = X^  and  X^(p-l)  = X2 ( t— l) . 

Proof : Let  X.^  = mXg.  construct,  as  in  the  previous  theorem, 
a cyclic  RM(t ,p,n,m,l)  design.  Take  Xg  copies  of  this 
design  to  obtain  a cyclic  balanced  RM(t,p,n,X1,X2)  design, 
as  desired. 

Even  though  shifts  prove  successful  in  constructing  bal- 
anced RM  designs  with  p < t,  this  is  not  the  case  for 
p = t.  Note  that  for  odd  t,  cyclic  balanced  RM(t,t,t,l,l) 
designs  do  not  exist.  This  is  so  because  in  this  case  we 
only  have  the  trivial  partition  = {1,2,3, .. .,t-lj.  Clearly, 

E j = = 0 (mod  t),  no  matter  how  we  order  the  elements 

J=1 

in  P1»  Nonexistence  is  now  assured  by  Proposition  3.2. 

The  construction  of  balanced  RMltjPjnjXi^Xg)  designs 
for  Xg  l X^  is  researched  presently.  Some  results  have 
been  obtained  and  they  will  be  presented  in  a subsequent  paper. 


4 . Estlmablllty  and  Estimation. 

A statistical  design  is  of  practical  use  if  at  least  the 
differences  between  various  parameters  of  interest  can  be  unbiasedly 
estimated  from  the  observations.  In  this  section  we  shall  be 
mainly  concerned  with  the  concept  of  estimability,  so  we  recall 
the  usual  definitions . 

Let  0 be  the  vector  of  parameters  and  Y be  the  vector 
of  observations . Suppose  we  assume  m effects  in  our  linear 
model.  Then  we  can  write: 

f>'  - 

where  0^  is  the  vector  of  parameters  associated  with  the 
ith  effect. 


Definition  4.1.  The  difference  between  two  parameters  is 
called  an  elementary  contrast.  A linear  combination  of  para- 


meters T,  CjYj  is  called  a (linear)  contrast  if 


T C j = 0 • 

J J 


Definition  4 .2 . A linear  function  of  parameters 
is  estimable  if  it  is  equal  to  the  expected  value  of  some 
linear  function  of  the  observations;  i.e.,  I* y is  estimable 
if  l#Y  = E(t'Y)  for  some  vector  t of  constants. 


Definition  4.3.  A design  is  connected  for  the  ith  effect 
if  all  the  linear  contrasts  of  parameters  of  the  ith  effect  are 
estimable.  A design  is  called  connected  if  it  is  connected 


1 

13. 

for  all  the  assumed  effects . 

We  should  mention  that  the  subspace  of  linear  contrasts 
of  t parameters  is  of  dimension  t-1  and  is  spanned  by  the 
elementary  contrasts . 

The  following  lemma  will  be  helpful  in  proving  the  next 
proposition. 

Lemma  4.1.  Let  S = {a^-a^,,  ^"^r+l*  a2~®r+2 * * * * ,at-l”ar+(t-l) ^ * 
where  ai  are  constants , and  let  (r,t)  = 1.  Then  ai-aj 
( i 4 J ) can  be  expressed  as  a sum  of  elements  of  S . (All 
the  indexing  is  modulo  t) . 

Proof:  Since  (r,t)  = 1,  the  Z-module  generated  by  r in 
Zt  equals  Zt  so  we  can  write  S = {aQ-ar,  ar”a2r»  a2r“a3r***** 
a(t-l)r-atr*  * 1 = kr  and  <3  = mr*  then  ai  - a^  = 

= ^akr'a(k+l)r^  + (a(k+l)r"a(k+2p+* * •+^a(m-l)r"amr^ * 

In  order  to  show  that  the  designs  constructed  in  the  pre- 
vious section  are  connected,  we  introduce  the  following: 

Proposition  4.1.  Assume  p < t . A cyclic  balanced 
RM(t,p,n,\1,X.p)  design  d is.  connected  for  direct  and  resi- 
dual effects  if  there  exist  two  pxt  subarrays  of  d such 
that  the  zeros  in  the  first  period  of  these  subarrays  are 
followed  by  a and  b respectively,  and  (b-a,t)  = ] 
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Proof;  The  necessary  conditions  n = X^t,  X-L(p-l)  = X2(t-l) 
together  with  p < t imply  X-^  > X2  > !•  The  design  d has, 
therefore,  at  least  two  pxt  subarrays.  Suppose  they  are  as 
follows ; 

(0  1 2 (t-1)  /0  1 2 (t-1) 

^a  a+1  a+2  ...  a+(t-l)  'b  b+1  b+2  ...  b+(t-l)  (*) 


where  a and  b satisfy  (b-a,t)  = 1.  Consider  the  entries 

(1,J),  ( 2 , j ) , (l,t+j),  (2,t+j)  in  (*).  We  have 

E(Yij)  = 4 + a±  + Sj  + TJ_1,  E(Ylft+Jj  = 4 + al  + ^t+J  + TJ-1’ 

E(Y2J)  = u + a2  + pj  + Ta+J_i  + Pj_i  and  £(Y2,t+J^  = 

= U + a2  + Pt+J  + + Pj_r  So 

E[?^Yij'Y2j  + Y2 , t+J_Yl, t+J^]  = Tb+J-1  “ Ta+j-l*  for  1 1 J 1 t* 

The  computations  in  the  indices  of  parameters  are  mod  t.  We 

have  estimated  t elementary  contrasts  for  the  direct  effects. 

The  previous  lemma  assures  us  that  all  the  elementary  contrasts 

for  the  direct  effects  are  estimable.  Let  T be  a linear 

rs 

A 

combination  of  the  measurements  such  that  E(Trg)  = Tr  “ Ts * 

Now  consider  the  entries:  (1,J),  (2,j),  (l,t+a-b+j), 
(2,t+a-b+j)  in  (*).  we  have  the  following  expectations; 

E<Y1J>  = M + ax  + Pj  + B(lf1,t+a.b+j)  - 

M + “1  + et+a-b+J  + Ta-b+J-l*  - » + ag  + Pj  + 

+ Ta+J-1  + PJ-1  and  E*Y2,t+a-b+J^  - M + a2  + Pt+a.b+J  + 
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+ Ta+J-1  + Pa-b+J-1*  Then 

E^Ylj‘Y2j+Y2,t+a-b+j‘Yl,t+a-b+J^  + ^a-b+j-1,  J-l1  = 

= Pa_b+ j_i"P j_]_ • This  holds  for  1 < J < t.  Again,  Lemma 
4.1  implies  that  there  exists  R^,  a linear  combination  of 
measurements,  such  that  ElR^)  = Pu  ” Pv* 

When  constructing  cyclic  balanced  RM  designs  it  is 
advisable,  for  purposes  of  estimation,  that  the  shifts  satisfy 
the  condition  (b-a,t)  = 1,  as  required  in  the  above  proposi- 
tion. 

Corollary  4.1.  The  designs  constructed  in  Theorem  3.1  are 
connected  for  direct  and  residual  effects . 

Proof : In  these  designs  (with  the  exception  of  X = 2 and 
p odd),  the  zeros  in  the  first  period  of  the  first  two  pxt 
subarrays  are  followed  in  the  order  of  application  by  l and 
2 respectively.  Proposition  4.1  assures,  therefore,  connected- 
ness (for  direct  and  residual  effects).  In  the  case  of  X-^  = 2 


and  p 

odd. 

the  difference  between  the 

same  two 

entries 

is 

(t-2). 

But 

(t-2 ,t ) =>  1,  since  in  this 

case  t 

is  odd . 

So 

Proposition  4.1  can  be  successfully  applied  to  this  case  as 
well.  Connectedness  for  the  residual  and  direct  effects  is 
therefore  established. 
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The  previous  proposition  assures  connectedness  for  direct 
and  residual  effects.  This  means  that  the  ranks  of  the  C- 
matrices  for  these  effects  is  at  least  t-1.  One  of  the 
results  that  follow  will  give  us  a little  more  insight 
into  what  these  ranks  can  be. 

Let  d be  a cyclic  balanced  RM(t ,p,n,X 1,X2 ) design 
with  design  matrix  X.  The  results  we  will  establish  in 
the  remainder  of  this  section  hold  even  when  the  design  d 
does  not  satisfy  part  (a)  of  Definition  2.1.  The  restriction 
p < t present  in  the  previous  section  won't  be  assumed  here. 
The  X'X  matrix  of  the  design  d is: 


T 

P 

a 

0 

M 

T 

VtxfVtXt 

XlJtxp 

P11  P12  ** 

P 

(p-lhjlt  t 

P 

P21  P22  •* 

* F2X1 

(p-UXji. 

X'X  = ct 

^l-^xp 

Jpxn 

Ul1pxl 

6 

Plnxn 

plnxl 

M 

pnlxl 

where  I is 

the  identity  matrix,  J is 

the  matrix 

with  all  the 

entries  equal 

to  1,  J*  is 

the  matrix 

with  first 

column 

entries 

all  zero  and  the  rest  of  the  entries  1.  The  symbol  laxi  stays 
for  the  column  vector  of  length  a with  all  the  entries  1.  The 
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IX  J* 


incidence  matrix  between  r and  p is  B1  * j^P.^  P12  •••  pi\ 
Similarly,  the  incidence  matrix  between  p and  p ia 

B2  = [Pgi  P22  P2X  1*  The  matrix  pij  (respectively  P2j) 
is  a txt  incidence  matrix  between  t (respectively  p)  and 

4*  h 

the  units  in  the  j subarray  of  the  design  d.  Since  each 
subarray  consists  of  shifts  of  (0,1,2, ... ,t-l ) it  will  be 
enough  to  find  the  first  row  of  P^,  the  other  rows  being 
cyclic  rotations  of  the  first.  There  are  exactly  p rows 
in  each  subarray  so  the  sum  of  the  entries  in  any  row  of 

will  be  p (all  1 < j < X-^).  Similarly,  the  sum  in  any  row 

of  P2J  will  be  p-1  (all  1 < j < X^.  If  we  let  P be 

the  following  txt  matrix: 


1 0 
1 0 

P = 

*0 
1 0 
1 0_ 

t“1  i 0 t 

then  Pi i = T a.P  where  P = P = I and  a. 's  are  non- 
J i=0  1 i 

t-1 

negative  integers  satisfying  E a,  = p.  Similarly 

1=0 

t-1  L 

P91  = E b.P  where  b.  are  nonnegative  inetgers  and 
0 i=0  1 

t-1 

E b,  = p-1.  If  the  design  d satisfies  condition  (a)  of 
i=0 

Definition  2.1,  which  is  the  case  with  all  the  designs  constructed 
in  Section  3,  the  integers  a^s  and  bi's  will  all  be  1. 

We  are  now  ready  to  introduce  a couple  of  useful  lemmas. 
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Lemma  4.2.  The  .joint  C- matrix  for  direct  and  residual 


effects  is: 


Ph1  - IjI/l/lJ  X2J'l2I  ‘ 


C(t>p)  = 


- - h T?o<K 


V-X2I  - 1 xl(p-i)i-  -^r-^  - PJ=V2/2J 


Proof:  Partition  the  X'X  matrix  as  follows: 


X1X1  X1X2 


X'X  = 


X2X1  X2X2 


, where  X'1X1 


P^lHxt  X2Jtxt  " Vtxt 


(P“l)x].Itxt 


lPll  P12  ***  PU1  ^iHxl 


XiX2  = 


I P21  P22  P2X ! 


nIpxp  Jpxn  nlpxll 


(p-l)X1ltxl 


X2X2  = 


plnxi  • The  matrix  C(t,p)  is 


obtained  as:  C(t,p)  = ” X1X2^X2X2^  X2X1  where 


It  is  a well-known  fact  that  the  Moore-Penrose  inverse 
of  a symmetric  matrix  is  symmetric,  commutes  with  the  matrix 
and  has  the  same  kernel  as  the  matrix.  The  following  lemma 
gives  a way  of  finding  the  Moore-Penrose  inverse  of  a poly- 
nomial matrix  in  the  txt  matrix  P introduced  earlier. 

Let  X be  an  indeterminate  and  f(X)  a (real  or  complex) 
polynomial  of  degree  at  most  (t-l)  in  P. 


Lemma  4.5.  The  Moore-Penrose  inverse  of  f(P)  = e a4P 
t_x  i=0  1 

is  f(P)~  = E b^P1,  where  b^^  = CQ+c-^w^+CgW2^. . .+ct_^w^ 

w is  a primitive  tth  root  of  unity  and  c^  equals  f(wi)“1 
if  f(w^)  4 0 and  0 otherwise. 

Proof:  The  txt  permutation  matrix  P satisfies  Pfc  = I 
and  has  in  fact  X^  - 1 as  minimal  polynomial.  Let's  in- 


troduce  the  following  useful  txt  matrix; 


20. 


4*  h 

where  w Is  a primitive  t root  of  unity.  The  matrix 
S is  unitary  and  it  is  easy  to  show  that  S_1PS  = D,  where 
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The  Moore-Penrose  inverse  of  f(D)  is 


f(D)~  = 


'y 


where  ci  equals  if  f(w*)  4 0 and  0 other- 

wise. Set  f(P)"  = Sf (d)”  S"1.  It  is  not  difficult  to 
check  that  f(P)~  defined  above  is  the  Moore-Penrose  in- 
verse of  f(P).  The  (i,J)th  entry  in  f(P)”  is 

1 ]r1  Vr  4 __ 

4 = t*  £ cvw  w . f(P)~  is  a polynomial  in  P,  because 

t j^3q  k 

tet1  = tr1  ck»kl  wkJ  . V ckwki  w«wk  wk  = ‘‘r1  ckwk<1+1)wk<J+1)  = 
k=0  k=0  k=0  K 

= te  j+1  (the  indices  in  e.^  are  computed  modulo  t). 

The  coefficients  b^  can  be  now  defined  as 

1 i ft— 1^1  - b— 1 . 

b.^  ® — (cq+c-jW  +...+ct_^w^  “ ' ) and  hence  f(P)“  =*  E b^P  . 


The  theorem  that  follows  reveals  some  properties  of  the 
C-raatrices  for  direct  and  residual  effects  obtained  from 
cyclic  balanced  RM  designs.  For  a positive  integer  d, 
let  ep(d)  denote  the  number  of  positive  integers  less  than 
d and  relatively  prime  to  it. 
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Theorem  4.1.  In  a cyclic  balanced  RM(t  jP^nA-^Xg)  design 

the  g-mat rices  for  direct  and  residual  effects  are  symmetric 

polynomials  in  p.  Their  ranks  equal  (t-1)  - E<p(d.)  £ex 

1 1 

some  (not  necessarily  same)  distinct  divisors  dj^  of  t. 

In  these  designs  contrasts  and  only  contrasts  of  direct  and 
residual  effects  may  be  estimable. 

Proof:  From  C(t,p)  in  Lemma  4.2,  we  obtain  at  once  the 
two  c-matrices: 

c(t)  = PXji  - i I1  - 

X1  J X1 
- <V  - V - 5 pijpb>  f(p,‘  <V  - ^ - 5 jfi  Vu’ 

\,(p-l)  , X1 

where  f(P)  = X -j^Cp-l) I ^ — J - | F P2JP^,  and 

J = 1 

# x x Xi(P-l)  i X1 

C(p)  - Xl(p-I)i-  -P5-;  - ? zx  - 
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Both  f and  g are  symmetric  polynomials  in  P with 
rational  coefficients.  In  virtue  of  Lemma  4.3,  we  can  set 
C(t)  = h(P)  and  C(p)  = k(P),  where  h(P)  and  k(P)  are 
themselves  symmetric  polynomials  in  P with  rational  coeffi- 
cients. The  dimension  of  the  kernel  of  h(P)  (respectively 
k(P))  is  equal  to  the  number  of  tth  roots  of  unity  which 
are  also  roots  of  h(X)  (respectively  k(X)).  Suppose 
w1,  a dth  root  of  unity  for  some  divisor  d of  t,  is 

J.  L_ 

also  a root  of  h.  Then  all  the  primitive  d roots  of 
unity  are  also  roots  of  h as  Images  of  w1  under  the 
Galois  group  of  automorphisms  of  Qtw1)  over  Q,  the  field 

±.  u 

of  rational  numbers.  There  are  qj(d)  primitive  d roots 
of  unity,  the  roots  of  the  cyclotomic  polynomial  of  w1 
over  Q.  We  now  show  that  1 is  always  a root  of  h.  As 
we  pointed  out  previously,  P2J  is  a polynomial  in  P whose 


t - 


entries  along  any  row  add  up  to  p-1.  Using  this  fact  in 

calculating  f(l),  one  gets  f ( 1)  = X,(p-l) - 1 

-L  tp 

- ^x(P-l)2  - X1(P-l)[l  - ~ - = 0.  The  same  is 


true  for  g,  the  argument  being  analogous.  Moreover, 
h(l)  = g(l)  - m(l)f(l)~  m#  (1)  = 0 since  both  g(l)  and 
(by  Lemma  4.3)  f(l)~  are  zero.  All  the  remarks  about  h 

hold  true  for  k as  well*  This  settles  the  statement  about 
the  ranks  of  the  C-matrices  in  the  theorem.  Having  always 
h(l)  =»  0 and  k(l)  = 0 means  that  the  row  sums  in  C(t) 
and  c(p ) are  always  zero.  This  implies  that  the  two 
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matrices  are  singular  and  hence,  of  rank  at  most  t-1.  We 
know  that  t#T  is  estimable  if  and  only  if  t'  is  in  the 
row  span  of  C(t).  But  t'  is  in  the  row  span  of  C ( r ) 
only  if  the  sum  of  the  entries  in  t is  zero  (i.e.,  only 
if  t'T  is  a linear  contrast).  The  same  holds  true  for  the 
residual  effects.  This  completes  the  proof. 

Let  the  four  txt  submatrices  of  c(r,p)  (see  Lemma 
4.2),  which  we  denote  by  A^j,  be  such  that 


C(t,p) 


A11  A12 


A12  A22 


Then 


are  polynomials  in  P.  By  evaluating  and 


A22  at  the  t roots  of  unity  it  can  be  readily  seen  that 
their  ranks  are  always  t-1.  By  deleting  the  last  row  and 
last  column  of  these  matrices  we  obtain  two  (t-l)  x (t-1) 
matrices  which  we  denote  by  Bi;l  and  B22<  Then 


are  generalized  inverses  of  A11  and  A22  respectively. 
The  usual  row  operations  in  the  normal  equations  lead  us  to 
the  following  result: 
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Proposition  4.2.  The  best  linear  unbiased  estimators  in  a 
cyclic  balanced  RM(t ^jnA^Xg)  design  for  the  estimable 
linear  contrasts  t'T  and  t' p are  respectively: 


t#T  = t'c(T)  [(x11-x#11x2(x#2x2)  x'2)-  a12a22(x^2-x5L2x2(x^x2)  xj)]y 

and 

t'p  = t'0(p)-[(xi2-X'12X11(X^X2)-Xj)-A21A-1(x'11-Xi1X2(X^X2)-x')]Y 

where  p = incidence  matrix  between  direct  effects  and 
units  in  the  jth  subarray. 

P2j  = incidence  matrix  between  residual  effects 
and  units  in  the  subarray. 


- PXXI  - p t PijPij 

0=1 

1 *1 

- V " V “ P ^ Plj^J 

. . x^p-D  , xi 

- - -Tp-J P2J^J 


C ( t ) “ A11  " A12A22A12 
c(p)  - A22  - AJ2A11A12 
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(xW  - 


-Lj 
pn  nxP 


pnJpxn  °pxl 

p^xn  °nxl 

°lxn  °lxl 


X1]L  = columns  in  the  design  matrix  corresponding 
to  the  direct  effects. 

X = columns  in  the  design  matrix  corresponding 
12 

to  the  residual  effects. 

Xg  = columns  in  the  design  matrix  not  in  X-^  or  X^>. 

Y =*  vector  of  observations  (ordered  the  same  way  as  in 
the  design  matrix) . 

Remark.  When  the  design  is  connected  C(t)“  and  C(p)~  can 
be  computed  the  same  way  as  A-^«  Otherwise  Lemma  4.3  can 


be  used 
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We  have  C(t)t  = TY  and  C(p)p  = RY  where  T and  R 
have  been  explicitly  given  above.  It  can  be  shown,  in  a 
straightforward  manner,  that  both  T and  R can  be  parti- 
tioned into  X-jP  submatrices  of  size  txt  each  and  each 
such  submatrix  is  a polynomial  in  the  matrix  P,  with  zero 
row  sums.  Let  these  partitions  be  T = £t^  . .. 

R " [R1  ^p]* 

Theorem  4.2 . In  a cyclic  balanced  and  connected  RM(t,p,n,X1,X2) 
design  if  the  b.l.u.e.  of  tq  " Ti  (respectively  pQ  - p^ 
is  Known,  then  the  b.l.u.e.  of  any  linear  contrast  of  direct 
(respectively  residual)  effects  can  be  derived  from  it . 


Proof  t Recall  that  for  a given  vector  l'  = (lltl2>  • • • >1^) 
its  shift  starting  at  i is  the  vector  (...  ltU1U2  • ••)'» 

i.  u 

denoted  by  l(t'),  in  which  is  in  the  i position. 

Indices  are  computed  modulo  t.  Suppose  l'C(r)r  = rQ  - t-^ 
where  l*  * . . .,<tt)  • Then  the  b.l.u.e.  of 

t0  - t ^ is  U'T^  ...  -C/ T^  p)Y.  Because  all  the  matrices 
Involved  are  polynomials  in  P,  - ^i+i  = 1(t#)c(T)t. 

Hence  the  b.l.u.e.  of  - Ti+1  is  (if-t'jT-^  ...  i(4' )T^p 
for  all  0 < 1 < t-1.  This  can  be  formulated  in  terms  of 
the  entries  of  the  design  as  follows:  The  first  row  of  treat- 
ments in  each  subarray  of  size  pxt  should  be  always  the 
same,  say  0,1,..., t-1.  Associate  the  weights  (coefficients) 
of  the  b.l.u.e.  of  - t1  with  the  corresponding  entries 


)Y 
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(cells)  in  this  design.  To  compute  the  b.l.u.e.  of 

- Ti+1,  leave  the  weights  associated  with  the  same  cells, 
but  do  cyclic  column  rotations  in  each  subarray  until  the 
cells  which  were  filled  by  0 will  be  filled  by  i.  The 
weight  associated  with  a certain  cell  will  then  become  the 
coefficient  of  the  observation  which  falls  in  that  cell  after 
these  cyclic  rotations  of  columns.  An  example  will  follow. 
Any  other  linear  contrast  of  direct  effects  can  be  expressed 
as  a linear  combination  of  the  t elementary  contrasts  con- 
sidered above.  The  b.l.u.e.  of  the  linear  contrast  is 
then  the  same  linear  combination  of  the  b.l.u.e.'s  of  these 
elementary  contrasts.  The  same  holds  true  for  residual 
effects . 

Let's  illustrate  by  a small  example  what  the  previous 
result  says.  Given  the  RM(3,2,6,2,l) 

0  12  0 12 

1  2 0 2 0 1 

b.l.u.e.  for  To  ” Ti 

^11  + *12  * 2*13  - *14  - *15  + 2*16  - 
- Y21  - Y22  + 2Y25  + Y24  + Y25  - 2Y26)  . 

b.l.u.e.  of  t2~t0  can  be  found  by  writing  the  design 


the 


The 

as 


■ 


2  C 1 
0 12 


2 0 1 
10  2 
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and  using  the  same  (unchanged)  weights  as  previously: 


T2  " T0  = 3\  13  + 

Y11  • 

2*12  ' 

Yl6  - 

y14  + 

2Y15 

i 

K"\ 

CM 

1 

*21  + 

2*22  + 

Y26  + 

y24  - 

2yE5 

The  b.l.u.e.  of  2tq  " ti  " t2  can  be  founci  by  subtracting 
the  b.l.u.e.  of  r2  - tq  from  that  of  rQ  - 

A list  of  partitions  for  values  of  t between  3 and 
20  is  attached.  From  each  partition  one  can  exhibit  the 
corresponding  balanced  connected  designs  as  shown  in  the 
proof  of  Proposition  3.2  and  the  two  examples  following  it. 

A comprehensive  bibliography  on  Repeated  Measurements  Designs, 
prior  to  1975,  can  be  found  in  A.  Hedayat  and  K.  Afsarinejad 
(1975). 
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Appendix 


Proof  of  Theorem  3.1. 

We  proceed  in  showing  that  the  conditions  stated  In 
Proposition  3.2  are  satisfied  for  all  the  partitions  listed 
after  the  statement  of  Theorem  3.1.  In  all  these  cases  it 
is  straightforward  to  check  that 
X1 

U P.  = (L,2,3>  • • . ,t-l>  (disjoint)  and  that  |P.  | =»  p-1 
i=l  1 1 

for  all  1 < i < 

k 

case  1.  (p  even).  To  check  the  condition  E s..  4 0 (mod  t) 

for  all  1 < J < k < (p-1),  where  P±  = ( (sitj ) < j < (p-i) 

is  the  P.^  defined  above,  we  note  that  sij  + s.^  j+1  = (t-l) 

for  all  odd  1 < J < (p-1)  and  all  1 < i < This  fact 

enables  us  to  classify  an  arbitrary  sura  i s..,  as  above, 

*=J 

in  one  of  the  following  four  patterns : 

(a)  [(i+JXx)  + (t-l)  - ( i+JXx) ]+...+  [ (i+kX-j^)  + (t-l)  - (1+kX^ 
for  all  1 < J < k < ^ ^ . If  for  some  j and  k 
this  sum  is  equal  to  mt  (for  some  me  z)  we  have 
(k-J+l)(t-l)  = rat.  Or,  (k-J+1)  =*  (k-J+l-m)t.  But 
2 < (k-J+l)  < < t,  which  leads  to  a contradiction. 
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(b)  [(i+JX-^  + (t-1)  - ( i+JX  ]+...+[  (i+fcX^  + (t-1)  - (i+kX1)] 

+ (i+(k+l)X1)  for  all  1 < j < K < . If  for  some 

J and  k this  sum  is  equal  to  rat,  we  get 
(-k+j-l+i+(K+l)x1)  =*  (ra-(k-j+l))t.  But  1 < i+j  = 

-k  + j-1  + i + (k+l)  < (-k+j-l+i+(k+l)x1)  = 

= (k+l)Xl  + i - 1 - (k-j)  < + Xx  = |X1  < (p-l)X1  = 

(t-l)  < t.  This  contradicts  the  fact  that 
(-k+J-l+i+(k+l)x 1 is  a multiple  of  t. 

(c)  (t-1)  - (i+(j-l)X1)  + [(i+jX1)  + (t-1)  - (1+JX-l ))  + 

...+  [(i+kX1)  + (t-1)  - (i+kX-j^)]  for  all 

0 < j-1  < k < ^--g  **■-.  This  sum  is  a multiple  of  t only 
if  (i+(  j -1  )X  ^K- J+2 ) is  a multiple  of  t.  But 

(i+( j-l)X1+k-J+2)  is  not  a multiple  of  t since 

2 < (i+(J-l)x1+k-J+2)  < \1  + -^^X-l  + + 2 < 

< l£^)Xl  + Ugl  < + Ugl  = (t-1)  < t. 

(d)  (t-1)  - (i+(J-l)Xl)  + [i+jX1+(t-l)-(i+j\1)]  + 

...+  [(i+KX1+(t-l)-(i+kX1)]  + i + (k+l)X;L  for  all 

0 < ( j-l)  < (k+l)  < ^-Pg2- « Proceeding  as  previously, 
this  sum  is  a multiple  of  t for  some  j and  k,  only 
if  (k-  j+2 ) (X  ^-l)  is  a multiple  of  t.  But 

2 < (k-J+2 ) (X  x-l)  < (k+l)Xl  < ^~X1  = t*'1)  < t. 

We  showed  that,  for  1 < i < X^,  the  conditions  on  P^'s 
required  in  Proposition  3*2  are  satisfied.  Since  i < X-^  has 


i 


not  been  used  in  any  of  the  previous  inequalities,  the  previous 

four  patterns  apply  also  to  P.  with  the  last  entry  removed. 

*1 

So  we  only  have  to  consider  the  cases  that  involve  the  last 

entry  of  P.  : i.e.,  t-1. 

*1 

(a1)  ( t — 1 ) + (t-l),  for  0 < l < — , where  4 is  the 

number  of  successive  pairs  of  entries  in  P that  add 

X1 

up  to  (t-l).  clearly  4(t-l)  + (t-l)  is  a multiple 
of  t only  if  (4+1)  is  a multiple  of  t.  But 

1 < 4+1  < < t . 

(b-^)  (t-l)  - 1 + 4(t-l)  + (t-l),  for  0 < 4 < and 

(j+4)  = ^ • Here  4 has  the  same  meaning  as  before. 

This  sum  is  a multiple  of  t only  if  (4+2+jX-j,)  is  also 

a multiple  of  t.  Since  2 < (4+2+jX^)  < (4+2)  + 

♦ - *K  s I ♦ s ? + s 

?t-l 

< — 2 — < t,  the  proof  for  p even  is  completed. 

case  2.1.  (p  odd  and  X^  > 2).  The  same  four  patterns 

arise  here  as  in  the  case  of  p even  and  they  are  ruled  out 

the  same  way.  Therefore,  we  only  have  to  consider  the  sums 

that  involve  the  last  entry  of  P.  . 

*1 

In  what  follows,  let  4 be  the  number  of  successive  pairs 

of  entries  in  P.  that  add  up  to  (t-l). 

X1 


I 


53. 

(a2)  4(t-l)  + + (t-l) ; for  all  0 < 4 < -^p-^2. 

Clearly,  4(t-l)  + + (t-l)  = mt  (for  some  m e z) 

only  if  (-4  + P2'1--X - l)  = (m-4-l)t.  But 

0 < i"1)  < ^~r~X  1 = < t# 

(b2)  (t-l)  - jX1  + t(t-l)  + + (t-l);  for  all 

0 < 4 < ^P2^ , and  j = ^Pg^  - 4 . If  this  sum  is  mt 
(for  some  m e z),  then  ( ^ Pg-1  - - j^X-^  “ 4-2  = (m-4-2)t. 
Since  -t  < -4  < - j)xi  - 4 - 2 < lP^ilx1  = i-t"1.)- 

< t,  we  can  only  have  f^Pg^-  “ ^)xi  = t “ 2 • Substitut- 
ing for  j = ^P-g^  - l,  we  get  \1  = But  is 

a positive  integer  if  and  only  if  4 = 0,  in  which  case 


case  2.2.  (p  odd  and  = 2).  P^^  is  defined  the  same  way 

as  in  the  previous  case,  so  we  should  only  examine  the  partial 
sums  in  P2 . But  the  partial  sums  in  Pg  with  the  first  entry 
removed  are  satisfactory  by  case  1.  Cases  involving  the  first 
entry  in  P2  should  only  be  considered.  There  are  two  possi- 
bilities: (a^)  (t-l)  + 4(t-l)  and  (b^)  (t-l)  + 4(t-l)  + 2J, 

where  0 < 4 < ^-Pg ?—  and  (in  (b^))  J = 4 + 1.  In  both  these 
cases  the  remainder  modulo  t is  4+1  and  1 < t + 1 < 


< t.  This  concludes  the  construction  of  the  designs  and  the 
proof  of  the  theorem. 
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t arbitrary,  p=2 : | (l) | (2) | (5) | ...  |(t-l)| 

t-5,  P-3 : | (1,3)  | (4,2)  | 

t-7,  P-3:  | (1,5) I (2,4) | (3,6) | 

t=7,  p=4:  |(1,5, 3)|(2, 4,6)| 

t=9,  p=3 : |(1, 7)|(2, 6) | (3,5) | (4, 8) | 

t=9,  P-5 : I (1,7, 3, 5) I (8, 2, 6, 4) | 

t=10,  p=4:  |(1, 8, 4)|(2, 7, 5)|(3, 6,9)| 

t=ll,  p=3:  |(1,9)|(2,8)|(3,7)|(4,6)|(5,10)| 

t-11,  P-6 : |(1,9,3,7,5)|(2,8,4,6,10)| 

t=13,  P-3 : |(1,11)|(2,10)|(3,9)|(4,8)|(5,7)|(6,12)| 

t=13,  p=4:  |(1,11,5)|(2,10,6)|(3,9,7)|(4,8,12)| 

t=13,  p=7 : |(1,H,3,9,5,7)|(12,2,10,4,8,6)| 

t=l5»  P-3:  |(1, 13)|(2, 12)|(3, 11)|(4, 10)|(5, 9)|6,8)|(7, 14)| 

t=15,  P-8:  |(1,13,3,11,5,9,7)|(2,12,4,10,6,8,14)| 

t-16,  p=4:  |(1,14,6)|(2,13,7)|(3,12,8)|(4,11,9)|(5,10,15)| 

t-16,  p=6:  |(1,14,4,11,7)|(2,13,5,10,8)|(3,12,6,9,15)| 

t-17,  P-3:  | (1, 15) | (2, 14) | (3, 13) | (4, 12) | (5, 11) | (6, 10) | (7, 9) i (8, 16) | 

t=17,  P-5:  |(1,15,5,H)|(2,14,6,10)|(3,13,7,9)|(4,12,8,16)| 

t-17,  P-9:  |(1,15,3,13,5,11,7,9)!(16,2,14,4,12,6,10,8)| 

t-19,  P-3:  |(1,17)|(2,16)|(3,15)!(4,14)|(5,13)|(6,12)|(7,11)| 

|(8,10)|(9,18)| 

t-19,  P-4:  |(1,17,7)|(2,16,8)|(3,15,9)|(4,14,10)|(5,3,11)|(6.12,18)| 

t-19,  P-7:  | (1,17,4,14,7,11) | (2,16,5,13,8,10) | (3,15,6,12,9,18) | 

t-19,  P-10:  |(1,17,3,15,5,13,7,11,9)|(2,16,4,14,6,12,8,10,18)| 
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